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ABSTRACT 

We studied the radio source associated with the ultraluminous X-ray source in 
NGC 5408 (Lx ~ 10^° erg s""'^). The radio spectrum is steep (index « —1), con- 
sistent with optically-thin synchrotron emission, not with flat-spectrum core emission. 
Its flux density (« 0.28 mJy at 4.8 GHz, at a distance of 4.8 Mpc) was the same in the 
March 2000 and December 2004 observations, suggesting steady emission rather than 
a transient outburst. However, it is orders of magnitude higher than expected from 
steady jets in stehar-mass microquasar. Based on its radio flux and spectral index, 
we suggest that the radio source is either an unusually bright supernova remnant, 
or, more hkely, a radio lobe powered by a jet from the black hole. Moreover, there is 
speculative evidence that the source is marginally resolved with a radius ~ 30 pc. A 
faint H II region of similar size appears to coincide with the radio and X-ray sources, 
but its ionization mechanism remains unclear. Using a self-similar sohition for the 
expansion of a jet-powered electron-positron plasma bubble, in the minimum-energy 
approximation, we show that the observed flux and (speculative) size are consistent 
with an average jet power « 7 x 10'^^ erg s^^ ~ 0.1ix ~ O.lLEdd, an age « lO'^ yr, 
a current velocity of expansion « 80 km s^^. We briefly discuss the importance of 
this source as a key to understand the balance between luminosity and jet power in 
accreting black holes. 

Key words: black hole physics - radio continuum: ISM — (ISM:) supernova remnants 
— X-ray: binaries — X-rays: individual: NGC 5408 X-l. 



1 INTRODUCTION 

Ultra-luminous X-ray sources (ULXs) are point-like, accret- 
ing X-ray sources with apparent isotropic luminosities > a 
few 10^* erg s~^, that is, greater than the Eddington limit 
of a stellar-mass black hole (BH) . Their natura remains un- 
explained. In particular, for the large majority of ULXs, two 
fundamental questions remain unsolved: whether the emis- 
sion is isotropic or beamed towards the observer (either as 
a mild geometrical beaming, or as a coUimated, relativistic 
jet); and whether the accreting compact objects are young, 
perhaps formed in recent starburst episodes, or are instead 
old relics of a primordial stellar population. 

Studying the ULX counterparts at other wavelengths 
may offer clues on their ages, mechanisms of formation, 
isotropic luminosities (through their effect on the interstellar 
medium) and sources of fuel. In particular, it has recently 
become clear that there is a strong conncction between X- 
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ray and radio properties of accreting BHs. This connection 
has been successfuUy applied to BH X-ray binaries in our 
Galaxy, providing a physical interpretation to their state 
transitions (Fender, Belloni & Gallo 2005). Similar studies 
for accreting BHs in nearby galaxies (e. g., Fender, Southwell 
& Tzioumis 1998) are, however, hampered by the faintness 
of the radio counterparts. 

The ULX in the nearby (d = 4.8 Mpc; Karachentsev 
et al. 2002) dwarf starburst galaxy NGC 5408 (henceforth, 
X-l) is an ideal test case for our study. Its position is known 
to < 0'.'5, from a series of Chandra observations. The source 
is among the brightest ULXs in the local Universe, with an 
isotropic X-ray luminosity « 1 x 10^" erg s~^ in the 0.3- 
12 keV band. Always detected in an active state also by 
Emstein, ROSAT and KMM-Newton, X-l has shown a pat- 
tern of flux variability over a few weeks or months (Kaaret et 
al. 2003), as well as on shorter timescales (Soria et al. 2004), 
inconsistent with an X-ray supernova remnant (SNR). Thus, 
there is no doubt that the X-ray emission is due to accre- 
tion onto a compact object. Most importantly, it is one of the 
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very few ULXs with a radio counterpart (Kaaret et al. 2003). 
The possibility that it is a background object (radio galaxy 
or quasar) can be ruled out due to the extremely high X-ray 
to optical flux ratio (Kaaret et al. 2003); a background blazar 
is also vcry unlikcly, bocausc their radio spectral indiccs 
are much flatter than determined for this source (Kaaret 
et al. 2003 and Section 2 of this work) . 

Two alternative scenaxios have been suggested for X-l 
(Kaaret et al. 2003; Soria et al. 2004). It could be a stellar- 
rnass microblazar (i. c, a microquasars with its rclativistic 
jet oriented along our line of sight). In this case, its ex- 
treme brightness would be due to Doppler boosting; its ra- 
dio emission would bc synchrotron from thc jet; its X-ray 
spcctrurn could bc produccd by synchrotron sclf-Cornpton 
or cxtcrnal Compton scattering of low-cncrgy photons (in- 
cluding UV/optical photons from the disk or the donor star) 
by the relativistic jet (Kording, Falcke & Markoff 2002; 
Georganopoulos, Aharonian & Kirk 2002; Gcorganopoulos 
& Kazanas 2003). Altcrnativcly, the X-ray emission could be 
isotropic, implying a mass > IOOM0 for the accreting BH, 
if the Eddington limit is not persistently violated. Apply- 
ing conventional phenomenological models, the 0.3-10 keV 
spcctrurn is dominatcd by a rapidly-varying powcr-law com- 
ponent with photon index F ~ 2.7, plus a non-varying soft 
thermal component or "soft excess" (fcT w 0.13 ke V), often 
found in ULXs. 

We re-observed the system in the radio, with a more ex- 
tended spectral coveragc (four instcad of two bands), to mca- 
sure its spectral index, and to determine whether there has 
been a long-term variation in the radio flux compared with 
the previous observations in Mareh 2000 (Stevens, Forbes & 
Norris 2002). Using these new results, we want to determine 
whcthcr thc radio emission is coming dircctly from the ac- 
creting BH (for example, from a relativistic jet) or is instead 
only spatially associated with it (for example, it may come 
from an undcrlying SNR). In cithcr casc, wc shall obtain 
important Information on thc gcomctry of emission or on 
the immcdiate environmcnt of thc ULX. More generally, we 
want to determine whether the associated radio and X-ray 
emissions in X-l are consistent with the correlations found in 
Galactic BH binaries, as a function of spectral state (Fender 
et al. 2004). 



2 ATCA RADIO STUDY 

Wc obscrvcd NGC5408 on 2004 December 9-11, with the 
Australia Telescopc Compact Array (ATCA), located at 
Narrabri, New South Wales. The array was in the 6-km con- 
figuration. Two sets of observations were carried out: simul- 
taneously at 1.344 GHz (20 cm) and 2.638 GHz (13 cm); and 
simultaneously at 4.840 GHz (6 cm) and 6.208 GHz (5 cm)\ 
with a continuum bandwidth of 128 MHz. Observations of 
the target galaxy were alternated with those of the primary 
flux calibrator PKSB1934— 638, and of the secondary phase 
calibrator PKS B1320-446. In total, wc had about 13 hr 
spent on source for each of the two configurations. 

The data were processed and analysed using the 

^ We chose to observe at 5 cm, rather than the more Standard 3 
cm, because of severe weather conditions, which made it impos- 
sible to preserve phase coherence at higher frequencies. 
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Figure 1. ATCA radio map at 4.8 GHz (6 cm). The X-ray posl- 
tion of the ULX is marked with a cross. The integrated 4.8-GHz 
flux from the ULX counterpart is R! 0.28 mjy. 



MIRIAD package (Sault, Tcubcn & Wright 1995). Images 
were cleaned with 1000 itcrations of gain 0.1. Thc clcarest 
image of the radio counterpart to the ULX is available at 4.8 
GHz (Figure 1), in which the source is clearly resolved from 
the nearby starburst and has a fiux density of 0.28 ± 0.07 
mJy (consistent with Kaaret ct al. 2003). The best-fit radio 
coordinatcs of thc source associated with the ULX, obtained 
from the peak of the emission in thc 4.8 GHz map, is R.A. 
= 14''03™■19^65, Dec = -41°22'58'.'7 with an error « 0'.'3. 
This is consistent (within 0'.'2) with the X-ray position of the 
ULX (R.A. = 14''03'"19".63, Dcc = -41°22'58'.'7; Kaaret et 
al. 2003) , which we also re-checked by analysing the Chandra 
observations available in the public archive. 

From the images at other wavelengths it was imme- 
diately apparent that the relative faintness of the target, 
combined with its proximity to the starburst region, was 
resulting in strong uncertainties in the measured flux den- 
sity. In order to quantify this, we have varied the beam size 
at each frcquoncy by adjusting Briggs' ROBUST parameter 
between -2 and +2, corresponding to a smooth transition 
from uniform (highest angular resolution, lowest signal-to- 
noise) to natural (lowest angular resolution, highest signal 
to noise). For each resultant map we have fit a point source 
model, with the coordinatcs fixed to thosc cstablishod from 
the 4.8 GHz map. The variation of the apparent source flux 
as a function of beam size, for each frequency, is presented 
in Figure 2. Beam size is clearly a strong effect at 1.3 and 
2.4 GHz, less so at 4.8 GHz and has no significant effect 
at 6.2 GHz: hencc, it is hard to confidcntly ascribe a well 
measured flux density to the radio counterpart of the ULX, 
particularly at low frequencies. A 4.8-GHz radio flux den- 
sity of (0.28 ± 0.07) mJy corresponds to a radio luminosity 
of (3.7 ±0.9) X 10^'' erg s~^. As noted above, the measured 
flux is consistent with the 4.8-GHz flux inferred from the 
March 2000 ATCA observations (Kaaret et al. 2003; Stevens 
et al. 2002), suggesting that the source has not varied sig- 
nificantly between the two epochs. 
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Figure 2. Changes in the apparent source flux as a function of 
beam size, in the four radio bands. At low frequencies, the source 
is contaminated by extended emission from nearby star-forming 
regions. See Section 2 for details. 



Between the highest two frequencies, the marginal de- 
tection at 6.2 GHz (in the naturally weighted maps, the 
point source fit at the position fixcd from the 4.8 GHz im- 
age is 0.12 ± 0.07 mjy) suggests at face value a very steep 
spectrum with spectral index an ^ —3^. This is ari cx- 
treme value for synchrotron emission: it may suggest that 
the break due to ageing of the electron populations is located 
between those two bands. We shall discuss — and esscntially 
rule out — this possibility in Section 3.6. Alternatively, it may 
be taken as a hint that we are starting to resolve the source 
at those higher frequencies, which would imply a character- 
istic size « 3" x 2", or « 50 x 70 pc at the assumed dis- 
tance of NGC 5408. At lowcr frcqucncics, contamination by 
the flux from the starburst starts to dominate. This makes 
determining the radio spectrum associated with the ULX 
problematic. 

In another approach to estimate the radio spectral in- 
dex at the ULX position and in the surrounding region, wc 
obtained radio images of the field by combining the 1.3 and 
2.4 GHz maps (Figure 3, top panel) and the 4.8 and 6.2 GHz 
maps (Figure 3, bottom panel). The radio source associated 
with the ULX clearly stands out as a discrete (and pos- 
sibly rcsolvcd) blob in the combined high-frcquency map, 
with hints of emission trailing back to the nearby bright 
starburst region. In the low-frequency map, there is clearly 
radio emission at the ULX location, apparcntly as the ond of 
an elongated region emanating from the starburst complex, 
but our spatial resolution is not sufficient for a more detailed 
analysis. We then took the ratio of the low-frequency and 
high-frequency maps (Figure 4). From this, we note that 
the bright starburst region is clearly strongly rcsolvcd out 
on its outer edges, which is why the spectral index in these 
regions appears steeper than —2. Instead, both the core of 
the starburst region and the site of the ULX are not resolved 
out, and have spectral indices w —1 between these two com- 



^ Throughout the paper, the spectral index is dcfincd such that 
F,^ ~ y" , and we use or ax to avoid confusion between the ra- 
dio and X-ray spectral slopes. Using the notation more commonly 
adopted in X-ray studies, the X-ray photon index F s — an, + 1. 




Figure 3. Top panel: radio map of the ficld obtained by combin- 
ing the 1.3 and 2.4 GHz data. The Chandra position of the ULX 
is marked with an error circle of radius C'e. Bottom panel: radio 
image obtained by combining the 4.8 and 6.2 GHz maps. 



bined bands. So, while a qualitative approach based on the 
mcasurod flux donsitics at four frequencics for the ULX is 
problematic, there is strong evidence that the radio emission 
at the site of the ULX has a steep spectral index, consistent 
with optically-thin synchrotron emission, and in agreement 
with the March 2000 observations (Kaaret et al. 2003). In- 
cidentally, tlic flux-ratio map also suggests that the radio 
emission in the core of the starburst region is dominated by 
SNRs (steep spectrum) rather than free-free emission from 
H II regions (flat spectrum). 
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Figure 4. Ratio of the high-frequency and low-frequency maps, 
showing that the ULX counterpart has a steep spectral index 
(«R ~ —1). The scale and orientation of the image is the samc 
as in Figure 3; the ULX position is marked with a cross. 



3 DISCUSSION 

3.1 Steady jet in a low/hard state? 

Given that the evidence for a spatially resolved radio source 
is still only marginal, we shall consider both the scenario of 
a point-like source (core emission) and that of an extended 
source. We start by considering the possibihty that the radio 
emission is directly produced by the accreting BH, on spa- 
tial scales comparable to the accretion disk. Radio emission 
from Galactic BH X-ray binaries is known to occur — with 
different physical properties — in two spectral states: the 
"low/hard" state (LHS) and the "very-high" state (VHS) 
(Fender et al. 2004). 

The LHS is characterised by a steady jet, with moder- 
ately low bulk Lorentz factor (7 < 1.4: Fender et al. 2004; 
Gallo, Fender & Pooley 2003). In a steady state, there is a 
"fundamental-plane" correlation between radio luminosity, 
X-ray luminosity and BH mass (Falcke, Kording & Markoff 
2004; see also Merloni, Heinz & DiMatteo 2003; Fender et 
al. 2004), of the form 



Lx CC L'SM^-^ 



where 



(17/12) - (2/3)qx 
(17/12) - (2/3)qr 



(1) 



(2) 



(Markoff et al. 2003). Taking this relation at face value — 
or, equivalently, using the empirical relation log Lr — 
0.60 log Lx + 0.78 log M + 7.33 (Merloni et al. 2003)— 
we would infer a BH mass ^ 10''Mq, well in excess of 
any estimates of other ULX masses, including the strong 
intermediate-mass BH candidate in M 82 (Strohmayer & 
Mushotzky 2003; Fiorito & Titarchuk 2004). 

However, there are at least two reasons why we rule out 
this result as unphysical. Firstly, the observed radio spec- 
tral index («r, « —1) is too steep to be consistent with 



an LHS steady jet, which is typically optically thick (core- 
dominated), with flat or slightly inverted radio spectrum 
(Fender et al. 2004). Secondly, the X-ray spectrum (Soria et 
al. 2004) is itself inconsistent with an LHS. The X-ray spec- 
trum in a LHS is dominated by a power-law of spectral index 
ax ~ —0.5 (i. e., photon index F « 1.5 as more commonly de- 
fined in X-ray literature). Instead, the observed X-ray spec- 
trum of X-l is much steeper (softer), with ax = —1.7 ± 0.2 
(i. e., r ~ 2.7). We conclude that the radio emission is not 
consistent with a steady jet from an accreting BH in the 
LHS. 



3.2 Flaring jet in the very high state? 

The X-ray spectrum (Soria et al. 2004) is consistent with 
the classical VHS of Galactic X-ray binaries, also called the 
"steep power-law" state, using the revised state classification 
of McChntock & Remillard (2003). Steady jets with a flat 
or inverted radio spectrum are sometimes observed in the 
"hard" sub-division of the VHS: this can be considered a 
high- luminosity extension of the LHS discussed (and ruled 
out) in Section 3.1. Stronger radio flares are observed in 
Galactic sources at the transition between the "hard" and 
"soft" subdivisions, as the sources become softer and cross 
the "jet line" (see Figure 7 of Fender et al. 2004). These 
flares are interpreted as discrete relativistic radio ejections, 
with a higher bulk Lorentz factor than in an LHS jet (7 ~ 
a few: Gallo et al. 2003). One explanation for the steep- 
spectrum radio emission is optically-thin synchrotron, from 
internal shocks in the jet, consequence of the sudden increase 
in the Lorentz factor during the state transition (Fender et 
al. 2004). 

Well-known examples of Galactic BH binaries with a ra- 
dio flaring behaviour in their VHS include GRS 1915-1-105, 
GRO J1655-40, and XTE J1550-564 (Table 1). If moved to 
the distance of NGC 5408, even at their peaks they would ap- 
pear two to three orders of magnitude fainter than the radio 
counterpart of X-l. Moreover, the radio to X-ray flux ratio 
is at least an order of magnitude higher in X-l, compared to 
those Galactic sources (Figure 5). Again, one might suggest 
that the enhanced radio luminosity, and radio-to-X-ray flux 
ratio are due to a much higher BH mass in X-l. 

However, even if X-l is interpreted (based on its X-ray 
properties) as a BH binary in the VHS, quantitative com- 
parisons with radio-flaring Galactic BHs are highly problem- 
atic. The VHS is a short-lived transient state, and a major 
ejection is a brief (rise time < 1 day), sporadic event (once 
every few months or years), as the source undergoes a hard 
to soft transition, associated with a change in the accretion 
disk inner radius. In our case, the radio flux was the same on 
the two occasions in which the source was observed, almost 
five years apart. Catching X-l during a major flare once 
would be a remarkable coincidence; observing it twice in ex- 
actly the same state would be too unlikely. Moreover, there 
is no evidence of any X-ray spectral state transitions over 
a series of Chandra and KMM-Newton observations span- 
ning the years 2001-2004. Hence, it is more likely that the 
observed radio emission is a persistent feature of this source. 

The only marginally viable alternative is that the source 
is persistently in the VHS, oscillating in a "dip-flare" cy- 
cle like GRS 1915-1-105, crossing the jet line back and forth 
in the process (Fender et al. 1999; Fender & Belloni 2004 
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Source 


Mbh 

(Mo) 


Date 
(MJD) 


Peak 5-GHz flux 
(mjy at 4.8 Mpc) 




Lyi/L-Edd 


Rcf. 


GRS 1915+105 (fir) 


14 


50750 


1.7 X 10^3 


0.6 


1.1 


F99 


GRS 1915+105 (osc) 


14 


50750 


2.6 X 10""' 


0.05 


1.1 


F99 


XTE J1748-288 


7 


50980 


1.7 X 10-3 


1.9 


0.1 


B06 


V4641 Sgr 


9 


51437 


1.2 X 10-3 


0.8 


4 


HOO, OOl 


GRO J1655-40 


7 


49580 


1.1 X 10-3 


1.0 


0.1 


HR95 


XTE J1550-564 


9 


51077 


2.0 X 10-"' 


0.3 


0.5 


W02 


XTE J1859+226 


7 


51467 


7.8 X 10-5 


0.2 


0.2 


B02 



Table 1. Parameters for the brightest radio outbursts in Galactic BH candidates. AU jet events listcd here occurred in the very-high- 
state, with both the jet power and the X-ray luminosity estimated to be ~ Eddington. Data from Fender et al. 2004; other references: 
F99 = Fender et al. 1999; B02 = Brocksopp et al. 2002; B06 = Brocksopp et al., in prep.; HOO = Hjellming et al. 2000; OOl = Orosz et 
al. 2001; HR95 = Hjellming & Rupen 1995; W02 = Wu et al. 2002. 
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Figure 5. Comparison between the X-ray and radio fluxes from 
NGC 5408 X-l and from some of the brightest Galactic micro- 
quasars, either in the LHS (steady jet) or in the VHS (flarings). 
In the on-line version of this plot, LHS datapoints are plotted in 
red, VHS in green. NGC 5408 X-l is a few orders of magnitude 
brighter, especially in the radio band. However, we argue that 
the comparison is not relevant, because the radio cmission from 
the X-l location comes either from an undcrlying SNR, or, much 
more likely, from a radio lobe. 



and references therein); this state is characterized by mod- 
erate X-ray and radio flaring on timescales of a few hundred 
seconds. Repeated flaring on similar timescales is indeed ob- 
served in the power-law component of the X-ray spectrum of 
X-l (Soria et al. 2004). Finding the source in that same dip- 
flare cycle in 2000 and 2004 (and in a similar X-ray spectral 
state consistent with the VHS in all observations between 
those two epochs) could perhaps suggest that this ULX has 
a different spectral-state behaviour than classical Galactic 
X-ray binaries. 

The peak radio flux density of GRS 1915+105 in the 
oscillatory state is three orders of magnitude lower than the 
radio flux of X-l, scaled to the same distance (Table 1). 
Taken at face value, the intensity of the radio emission from 
X-l in 2000 and 2004 make this ULX very peculiar, in com- 
parison with typical Galactic BH X-ray binary, in which the 
radio emission is believed to come from a compact jet (Fig- 
ure 5). However, we must also consider the possibility that 
the radio source associated to X-l is enhanced by Doppler 
boosting, as we shall discuss in Section 3.3. 



3.3 Relativistically-beamed emission? 

Relativistic Doppler boosting is a possible mechanism to rec- 
oncile the brightness of X-l with those of Galactic BHs, as 
originally suggested by Kaaret et al. (2003). The amplifi- 
cation of the rest-frame emission due to Doppler boosting 
~ S''-"^, where S = 7-^(1 - /3cos6')-^ and 2 «S fc 3 de- 
pending on whether it is a continuous jet or sporadic ejec- 
tions. Having already ruled out (Section 3.2) both a steady 
jet in the LHS and a single major flare, we only need to esti- 
mate the Doppler boosting required for a persistent dip-flare 
cycle. In GRS 1915+105, the Lorentz factor 7 in this state 
is > 2 and perhaps as high as 5 (Fender et al. 2004); at a 
viewing angle »; 65° , the observed 5-GHz flux is 50 mJy 
at 11 kpc, corresponding to ~ 0.3 fiJy at 4.8 Mpc. To en- 
hance this value by a factor ~ 10^, a viewing angle 9 < 20° 
(for fc = 3) or < 15° (for A; = 2) is required, for a source 
identical to GRS 1915+105 and 7 = 5. AUowing for an in- 
trinsic radio luminosity an order of magnitude stronger than 
in GRS 1915+105 (justified by the fact that the apparent 
X-ray luminosity is also a factor of 10 higher), only 9 < 30° 
(for fc = 3) or ^ < 25° (for k = 2) are required. As an 
aside, we note that in order to allow for persistent dip-flare 
ejections with Lorentz factors « 5, the microquasar would 
have to reside in a large low-density bubble, otherwise the 
jet-ISM interaction would rapidly reduce the jet speed. 

For the X-ray emission, strong beaming is disfavoured 
by the presence of a soft, non-variable thermal component 
at fcTbb = 0.13 ± 0.01 keV, with unabsorbed luminosity 
3 X 1039 erg 3-^ in the 0.3-1 keV band (Soria et al. 2004). 
It is still unclear whether this soft component, typical of 
many ULXs (MiUer, Fabian & Miller 2004; Feng & Kaaret 
2005), comes from the accretion disk or from a downscatter- 
ing corona or outflow (e. g., Laming & Titarchuk 2004; King 
& Pounds 2003), or perhaps from a combination of reflection 
and absorption (Chevallier et al. 2006; Crummy et al. 2006); 
but in any case, it is difflcult to reconcile with relativistically 
beamed emission. Thus, it is more likely that at least the soft 
component of the X-ray emission is isotropic, requiring a BH 
with a mass > IOOM0 to satisfy the Eddington limit. It is 
possible, instead, that the power-law component is beamed, 
for example produced by inverse-Compton scattering of disk 
photons by the jet (e. g., Bosch-Ramon, Romero & Paredes 
2005; Georganopoulos et al. 2002). 

We conclude that Doppler boosting by a relativistic jet 
can in principle explain the observed radio flux and is not 
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inconsistent with the X-ray properties. However, we consider 
this scenario somewhat contrived because it relies on the 
source spending most of its time in a "persistent" VHS with 
steep-spectrum dip-flare osciUations over many years. This 
is unlike the typical spectral state behaviour of Galactic X- 
ray binaries and microquasars; no other sources with the 
same behaviour have been found yet. The other difhculty is 
that we have at least speculative evidence to believe that 
the radio source is extended over > 30 pc (Section 2), which 
would rule out point-like core emission. 



3.4 Supernova remnant? 

In all the scenarios discussed — and in most cases ruled out — 
so far, the radio emission would be coming from the inner- 
most region around the BH, and would be directly related 
to its current state of accretion. We shall now consider the 
possibility that the emission comes, instead, from a larger 
region in the interstellar medium around the BH. 

The association with free-free radio emission from 
a star-forming region was already ruled out (Kaaret et 
al. 2003) due to the observed steep spectrum, inconsistent 
with thermal bremsstrahlung. A steep optically-thin syn- 
chrotron spectrum would instead be consistent with an un- 
derlying shell-like SNR: remnants in this class have typi- 
cal radio spectral indices — 1 < or < —0.4. On the other 
hand, we can rule out fiUed-center (plerionic, or pulsar-wind- 
nebula) SNRs — such as the Crab — because their spectral 
indices are much flatter (—0.3 < or < —0.1; e. g., Weiler 
& Panagia 1978). Our observed index or ~ — 1 would be 
among the steepest, but not unique. For example, at least 
three, and possibly five of the 14 shell-like SNR candidates in 
NGC300 have a spectral index < -1 (Pannuti et al. 2000); 
three of the 15 shell-hke SNR candidates in NGC 6946 have 
an index between —1.4 and —1.0 (Hyman et al. 2000); two 
of the 53 SNRs in M 33 (including both plerion and shell-like 
SNRs) have an index ^ -1.1 (Duric et al. 1995). 

The radio flux of our source in NGC 5408 0.28 mjy 
at 4.8 GHz) would also be among the highest for SNRs in 
nearby galaxies, but not unique. The brightest radio SNRs 
in NGC 300 would only reach a flux density ~ 0.06 mJy at 
the distance of NGC 5408. However, 6 of the 15 radio SNR 
candidates in NGC 6946 would have a 5 GHz flux > 0.2 mJy 
(Hyman et al. 2000); the brightest SNR in that galaxy would 
have a flux ~ 0.5 mJy. All five of the candidate SNRs in 
NGC 4258 (Hyman et al. 2001) are also brighter, with fluxes 
« 0.6-0.7 mJy at the distance of X-l. In NGC 7793, four of 
the five brightest radio SNRs would have fiuxes ~ 0.08-0.12 
mJy at the distance of NGC 5408 (Pannuti et al. 2002); how- 
ever, the brightest candidate SNR in that galaxy, NGC 7793- 
S26, has a spectral index qr — —0.9 ±0.2 and would have a 
5 GHz flux « 0.6 mJy. The brightest shell-like SNR in our 
Galaxy, Cas A, would be seen with a 5 GHz fiux of 0.3 
mJy (Reynoso & Goss 2002) and a diameter of « 0'.'2. 

One explanation for such exceptionally bright SNRs is 
that they are produced by very massive stars with strong 
winds; hence, the explosion occurred in a dense circumstellar 
environment (Dunne, Gruendl & Chu 2000). Alternatively, 
the remnants may be the result of multiple SN explosions, 
perhaps inside a star cluster (Pannuti et al. 2002). Either 
scenario would be consistent with a direct physical asso- 
ciation between the SNR and the process of formation of 



Figure 6. ATCA radio flux density contours at 4.8 GHz, over- 
plotted on a B image from Subaru (top panel) and an Ha image 
from the ESO 3.6m telescope (bottom panel), on the same spatial 
scale. 

the ULX X-l. In addition, a relatively young age may play 
a role. In summary, a bright shell-like SNR would be con- 
sistent with the measured fluxes and spectral index, would 
explain why the flux has been approximately the same be- 
tween 2000 and 2004, and would be a natural candidate for 
the formation of the ULX. If, as we speculate, the radio 
source in NGC 5408 is marginally resolved with a diameter 
~ 50-70 pc, we can apply the purely empirical correlation 
between age and diameter, found in Galactic SNRs (Xu, 
Zhang & Han 2005): we obtain an age > 10* yr, with a 
best-fit age ~ 10^ yr (but with a large scatter). Thus, the 
source would be too large to be a very young or histori- 
cal SNR (like Cas A). In Sections 3.5 and 3.6 we shall try 
to find more stringent tests to the SNR scenario, based on 
the minimum-energy condition, and suggest an alternative 
explanation. 

3.5 Hq emission: SNR or X-ray photoionization? 

The SNR scenario was initially rejected (Kaaret et al. 2003) 
because there appeared to be no evidence of optical emission 
lines at the X-l position (upper limit of ~ 7 x 10~^^ erg s~^ 
cm^'^ derived for Ha by Kaaret et al. 2003, from a broad- 
band HST/WFPC2 image in the F606W filter). Ahhough 
the ULX appears to be well separated from the starburst re- 
gion, in broad-band optical images (e. g., the _B-band Subaru 
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Figure 7. Close-up view of the region around X-l, in the 
continuum-subtracted Ha image from the ESO 3.6-m telescope. 
A circle (green in the on-hne version) marks the Chandra posi- 
tion; its l"radius combines the rms error in the astrometric so- 
lution for the opticai image with the astrometric uncertainty in 
the Chandra position. ATCA 4.8-GHz radio flux density contours 
are overplotted (cyan, in the on-line version) : the contours are for 
fluKCS of 0.054, 0.081, 0.12 and 0.18 mjy beam"!. 



image in Figure 6, top panel; and the HST/WFPC2 image 
in Kaaret et al. 2003), a narrow-band Ha image (Figure 
6, bottom panel) shows instead that most of the galaxy is 
filled by a tangled web of HII filaments and knots. For this 
study, we downloaded and analysed public-archive Ha and 
R-continuum images taken from the ESO 3.6m telescope 
(EFOSC camera) by M. PakuU and L. Mirioni. 

In particular, an Ha-emitting knot is visible at the ULX 
position (Figure 7) (as noted in Mirioni 2003, PakuU & Miri- 
oni 2003): it has a slightly elongated shape, with a size of 
« 50 X 70 pc, matching the upper limit (and possibly the 
true size) of the radio-emitting blob. We checked that the 
slight displacement between the Ha knot and the radio and 
X-ray positions is within the error in the astrometry of the 
opticai image 1"); therefore, the Ha, radio and X-ray 
emission are consistent with being coincident. The Ha knot 
is connected to the main starburst region through a net- 
work of filaments and cavities (Figure 7). We estimate a 
flux fua = (5 ± 2) X 10"^'' erg cm"^ s~^ from that knot, 
corresponding to a luminosity Lua ~ (1-5 ± 0.4) x 10^^ erg 
s~^ in the line (or ~ 5 x 10*^ Ha photons s~^), that is 
~ 10"^ times the total Ha emission from NGC 5408. The 
peak surface brightness Sua ~ 1.5 x 10^^^ erg cm~^ s~^ 
arcsec"'^. 

If the Ha line emission is due to recombination, its es- 
timated luminosity requires a flux of ~ 1-2 x 10''* ioniz- 
ing photons s~^ intercepted by the nebula (Osterbrok 1989; 
cf. also similar calculations for other photo-ionized nebulae 
in Pakull & Angebault 1986, and Soria et al. 2005). This is a 



relatively modest requirement. Continuum Subaru and HST 
images show three massive stars near the peak of the radio 
and Ha emission, with absolute B magnitudes « — 6 and 
bluish colors, B — V<0 (Kaaret et al. 2003; Copperwheat et 
al. 2006, in prep.). They are consistent with main-sequence 
O stars, or, more likely, blue supergiants. Three such stars 
with initial masses ~ 2OM0 would provide enough ionizing 
photons to explain the Ha emission (Sternberg, Hoffmann & 
Pauldrach 2003; Schaerer & de Koter 1997). As a consistency 
check, the size of a Stromgren sphere at a characteristic tem- 
perature T « 10* K is 7?s ~ 30 G^i^n^^^^ pc, where Q48 is 
the photo- ionizing flux in units of 10** photons and tih 
is the hydrogen number density in units of 1 cm~^. This is 
in agreement with the observed size of the H II region at the 
ULX position'^. 

We stress that this Ha nebula is not at all unique or 
remarkable for this galaxy, and is fainter than many other 
knots and filaments (Figure 6), which are probably due to 
OB photo-ionization or local density enhancements in the 
outflowing warm gas. It is only its association with strong 
X-ray and synchrotron radio emission that makes it remark- 
able. Hence, we also need to consider the possibility that 
the local enhancement in the Ha emission is directly caused 
by the ULX. We used XSPEC to re-analyze the XMM- 
Newton X-ray spectra of the ULX from 2001-2003 (Soria 
et al. 2004), using variable-abundance, ionized absorbers in 
addition to the Galactic line-of-sight term. We found that 
intrinsic absorption (column density ~ 1 x 10^^ cm~^) re- 
moves ~ 1.7 X 10** photons s~^ with energies above 13.6 eV, 
assuming a thermal Comptonization model (bmc) for the X- 
ray spectrum*. Plausible local interstellar gas densities nu ~ 
a few cm~"^ are consistent with a significant fraction of the 
intrinsic absorption coming from < 30 pc from the ULX. 
Hence, photoionization from the ULX may contribute a sim- 
ilar amount of H-ionizing photons as the few OB stars in the 
surrounding region. 

Finally, part of the Ha emission may be due to shock- 
ionization from an SNR. Most identified radio SNRs in 
nearby galaxies are associated with much brighter Ha emis- 
sion, from ~ a few 10~^^ erg s~' cm~^ up to ~ a few 10~^* 
erg s~^ cm~^ at the distance of X-l (Hyman et al. 2000; Pan- 
nuti et al. 2000; Deharveng et al. 1988). This is partly due to 
a selection bias (e. g., Pannuti et al. 2000): unresolved, non- 
thermal (steep-spectrum) radio sources without opticai line 
emission would generally be discarded as background AGN. 
However, radio synchrotron and Ha emission in SNRs come 
from different processes, and are not strongly correlated; in 
fact, there are a few radio-bright candidate SNR in nearby 
galaxies, with Ha emission below the level inferred for the 
source in NGC 5408: for example, some sources in the Large 



^ Note that the thermal (free-free) radio emission associated with 
a photoionized H II region of such size and luminosity is 500 
times fainter than the non-thermal emission measured at the ULX 
position (Dickel, Harten &; Guli 1983; Lazendic, Dickel &; Jones 
2003). Hence, the thermal contribution is absolutely negligible in 
our case, and does not afTect our estimate of the radio spectral 
index. 

* Power-law plus blackbody models also provide good fits to the 
X-ray spectrum, but are less suitable for an extrapolation to the 
UV region, since they give unphysically high fluxes, unless the 
power-law component is truncated at low energies. 
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MagcUanic Cloud (LMC) (Figure la in Filipovic et al. 1998) 
^. Othcr examples of Ha-faint SNRs are found in NGC300 
(Table 4 in Payne et al. 2004). 

In summary, we conclude that the SNR scenario for the 
radio synchrotron cmission rcmains viablo, although prob- 
ably not the most likely. We shall also see in Section 3.6 
that the radio source would require am explosion energy un- 
usually high for a normal SN. In the optical band, more 
narrow-band imaging and spectral studies are needed to dis- 
tinguish bctwcon the main compcting sccriarios (H II region, 
X-ray ionized nebula, shock-ionized SNR). Specifically, the 
most important constraints will come from the presence or 
absence of nebular [O III] and Hc II emission, and from the 
[S II] /Ha and [Ol]/Ha line ratios. 



3.6 Radio lobes? 

Finally, we consider the possibility that the optically-thin 
synchrotron radio emission comes from extended radio lobes 
associated with the ULX. While core emission from a ra- 
dio jet (Sections 3.1-3.3) traces the instantaneous accretion 
power and spectral statc of the BH, emission from radio 
lobes dcpends on the accretion activity integrated over the 
source lifetime (e. g., Scheuer 1974), even if the core were not 
currently active. Hence, we do no expect significant variabil- 
ity over a few years, in agreement with the observations. 

Almost all radio-loud AGN show a lobc structure, which 
dominates the radio emission in FR II sources. On the other 
hand, while most X-ray binaries are thought to release part 
of their accretion power in a jet, only a few sources have 
been found so far with persistent lobe-like emission: Cyg X-l 
(Gallo ot al 2005); Gir X-l (Tudosc et al. 2006); the Galactic 
microquasars lE 1740.7-2942 and GRS 1758-258 (Mirabel 
et al. 1992; Hardcastle 2005). In addition, in SS433, a 
radio-lobc structure may cocxist with an undcrlying SNR 
(Downes, Pauls & Salter 1996; Dubncr et al. 1998); it has 
been noted that SS433 may be classified as a ULX if we 
could see it face-on (Fabrika & Mescheryakov 2001). The 
size of the steep-spectrum radio-lobe structure in SS433 
is « (50 X 100)(d/3 kpc) pc, very similar to the source in 
NGG 5408; however, the radio luminosity density is an order 
of magnitude lower. 

The main difference between radio galaxy and micro- 
quasar lobes is due to the environment: while some quanti- 
ties, such as all charactcristic timc-scales, length-scalcs and 
energics, scale as Mbh, others, such as the density of the am- 
bient medium, do not (e. g., Heinz 2002). In particular, radio 
lobes in stellar-mass microquasars tend to be much more 
over-pressured with respect to the external medium; there- 
fore, thcy have a shorter lifetime due to adiabatic expansion 
(Hardcastle 2005; Hardcastle & Worrall 2000). Nonetheless, 
Palcull & Mirioni (2002) discovered that many ULXs are in- 
deed surrounded by laxge optical nebulae possibly powered 
by the interaction of a jet with the interstellar medium. 

The spectral index and total radio power alone do not 



^ For comparison with the plot in Filipovic et al. 1998, the Ha 
flux and 4.8 GHz radio flux density from the source in NGC 5408 
would be fs 5 X 10~^^ erg cm~^ and w 2.5 Jy, respectively, 
translated to the distance of the Large Magellanic Cloud. 



give us enough Information to determine the physical pa- 
rameters of the lobes — in particular, the energy density of 
the particles and the magnetic field. However, if the size of 
the lobes (hence, their surface brightness) is known, one can 
minimize the total energy density to get a rough constraint 
on the particle and field energy and pressure (e. g., Pachol- 
czyk 1970). In our case, the radio source may be marginally 
resolved (Section 2), although its detailed spatial geometry 
is unknown. For simplicity, we take a spherical lobe with a 
radius of 30 pc (~ l'.'3), a power- law distribution of rela- 
tivistic electrons with spectral index —3 (corresponding to 
an = —1), and we assume that the electron distribution 
cxtends between Lorcntz factors of 1 and 10* ^. We also as- 
sume a magnetic field randomly distributcd with respect to 
the line of sight, and a volume fiUing factor « 1. See Ap- 
pendix A for a brief summary of the formulae we used to 
estimate the minimum-energy magnetic field. 

The matter content of a quasar or microquasar jet is 
still a hotly debated issue (Romero 2005). For example, the 
Galactic microquasar LS 5039 has been variously modelled 
with leptonic (Bosch-Ramon et al. 2005; Dermer & Bottcher 
2006) or hadronic jets (Romero et al. 2003; Aharonian et 
al. 2005). An clectron-positron plasma was suggestcd for 
the jet in Nova Muscae (GRS 1124-684) (Kaiser & Han- 
nikainen 2002). On the other hand, a significant hadronic 
contribution could explain why synchrotron emission is not 
detected within the free-free emitting bubble of Cyg X-l 
(Gallo et a. 2005). Hadronic jet models are also favoured for 
the microquasars SS 433 (MigUari, Fender & Mendez 2002), 
and LS I -1-61303 (Romero, Christiansen & Orellana 2005). 
A leptonic model was instead suggested for the powerful 
quasar 3C 279 (Wardle et al. 1998), based on the detection 
of circular polarization. 

From the observed 4.8 GHz flux, for an electron- 
positron plasma, we obtain: a minimum energy density 
« 9.4 X 10"^" erg cm^"^ for a random magnetic field B ~ 0.11 
niGauss^; a mininmm pressure Pin ~ 5.4 X 10^^'^ dynes 
cm~^; and a total energy ~ 3.3 x 10"'^ erg. Alternatively, 
taking two lobes with radii of 23 pc = 1", the minimum 
pressure is « 6.2 x lO"'^" dynes cm~^ and the total energy 
is « 2 X 1.7 X 10^^ erg. Most quantities dcpend strongly on 
the adopted value of or. Taking or = —0.7 as a strong up- 
per limit, we infer a total energy > 3.5 x 10^" erg for a single 
30-pc electron-positron bubble, or > 3.6 x 10^° erg for two 
23-pc lobes. 

If energy and pressure in the synchrotron-emitting bub- 
ble are dominated by relativistic protons, the total energy 
can be up to a factor of 10 higher than in the leptonic case, 
for the same sets of parameters. For example, if hadrons 



^ For a random magnetic field Ri 0.1 mGauss, the synchrotron 
cooling timescale of relativistic electrons is RJ 10^ X (lO'^/'y) yr, 
where 7 is the Lorentz factor. Talcing instead an upper limit 7 ~ 
10^ for the electron distribution does not significantly affect our 
results: for our rathcr stccp spectral indcx Ri —1, the minimum 
energy dcpends more strongly on the low-cnergy cut-off of the 
electron distribution. For cxamplc, a lowcr cut-off at 7 ~ 10 would 
rcduce the estimated energies by a factor of 3. 
^ For the adopted spectral index or = — 1, the minimum-energy 
condition also implies energy equipartition, while for a spectral 
index = —0.7 the energy density of the magnetic field would be 
Ri 85% of the particle energy density. 



A ULX microguasar in NGC5408? 9 



carry ~ 99% of the energy*, the minimum- oncrgy condition 
implies a magnetic field B « 0.35 mGauss and a total en- 
ergy « 3.3 x 10^^ erg for a 30-pc synchrotron bubble with 
a spectral index «r = — 1; our limit «r < —0.7 implies an 
energy > 4.3 x 10*^ erg. Honcc, the leptonic-jet assumption 
can be used as a conservative lower limit to the energy in 
the bubble. 

On the other hand, relativistic hadrons do carry most of 
the energy in SNRs (known sources of cosmic rays). There- 
fore, if the radio source in NGC 5408 is interpreted as an 
SNR, we need to use the higher range of energy estimates. 
Energies ~ 10®^ erg are implausibly high for a single "Stan- 
dard" SN. FoUowing this argumont loads us to concludc that 
a SNR interpretation is less likely than a radio lobe scenario, 
but we cannot rule out the possibility of a more energetic 
hypernova remnant. Conversely, we could assume that the 
source is an SNR with a Standard energy w 10®^ erg, and in- 
fcr its sizo from the minimum-oncrgy condition. It turns out 
that the remnant should bc ~ 10-15 pc in radius, for the ob- 
served flux and spectral indcx. Hence, deeper high-frequency 
radio observations will be crucial to unequivocally determine 
the size of the radio source (or at least a more stringent up- 
per limit) and tost the SNR scenario. 

One might suggest that the injection spectral index was 
flatter (qr w —0.5), and the steeper index currently ob- 
served at high frequencies may simply be due to an ageing 
population of relativistic electrons inside the expanding bub- 
ble ( "spectral ageing" ) , as the higher-energy olcctrou popu- 
lation is depleted by synchrotron losses. If this is the case, 
the total energy in a 30-pc bubble would be lower, w 1 x 10®^ 
erg (corresponding to a magnetic field B « 0.07 mGauss), 
consistent with a normal SN. The presence of a spectral 
break in the radio continuum would, in principle, provide 
another criterion to distinguish between an SNR and a jet- 
inflated lobe. For an SNR, we expect an exponential cut-off 
in the spectrum abovc the break frequoncy (Jaffc &: Per- 
ola 1973); for a jet bubble, charactcrized by continuous in- 
jection, the spectrum would bc a brokcri power law with an 
index «r = q;r — 0.5 above the break frequency (Kardashev 
1962). 

The problem with the spectral break scenario is that 
a relatively old age is required for the break to occur 
at frequencies as low as ~ 3-4 GHz. FoUowing Murgia 
(2003), and using the previous order-of-magnitude esti- 
mate of the magnetic field, we obtain a characteristic age 
^8x10^ {B/0.1 mG)"^-'''(i^br/4 OMz)""-"' yr for an cxparid- 
ing 30-pc bubble. SNRs of that age are unlikely to be still 
so bright. This age is also longer than the characteristic age 
we shall estimate for the radio- lobe scenario (Section 4). 
However, even in this case we cannot entirely rule the break 
scenario, for example if B is well above equipartition in a 
magnetically dominated outflow. 



4 SELF-SIMILAR SOLUTION FOR RADIO 
LOBES 

In Section 3.6 we have treated the size of the bubble sim- 
ply as an arbitrary input parameter. We can learn more 
about the time evolution of a jet-inflated bubble if we ex- 
press its radius sclf-consistcntly as a function of the total 
energy injected into it, thus coupling a radial oxpansion 
equation with the energy density equation. To do so, we 
describe the supersonic (but sub-relativistic) expansion and 
radio emission from the bubble with a set of self-similar so- 
lutions (Castor, McCray & Weaver 1975; Fallc 1991), of- 
ten adopted to model radio lobes in AGN and Galactic 
micro-quasars (Heinz, Reynolds & Begelman 1998; Willott 
et al. 1999; Jarvis et al. 2001; Heinz 2002). We assume (Heinz 
et al. 1998): a spherical, uniform radio cocoon of radius r e', a 
constant density of the external ISM (an approximation suit- 
able to microquasar environments); relativistic gas inside the 
radio bubble (adiabatic index = 4/3); non-relativistic gas in 
the shell outside the contact discontinuity (adiabatic index 
— 5/3) ; equipartition of field and particle energies in the co- 
coon; non-radiative shocks; a radio spectral index or = — 1. 
For simplicity, we use a density p = 1.7 x 10~^'^nH, valid 
for a pure hydrogen gas; in any case, the main results from 
the self-similar solutions do not have a strong dependence 
on the mean atomic weight. 

The self-similar equations governing the expansion of 
the bubble are (c. g., Heinz et al. 1998): 

rc ~ 4.02 X 10*n-^/^£^/«^^/^ (3) 

Ct = Etot, (4) 

where Etot is the total energy injected by the jet into the 
radio cocoon, and £ is the time-averaged jet power, assum- 
ing that it is either approximately constant, or with a duty 
cycle much shorter than the time t. The emitted flux density 
at a frequency v is then 




See Appendix B for an outline of the derivation of Eq.(5). 
We foUowed Willott et al. (1999) to obtain the normal- 

ization of the rclations between jet power and age and radio 
flux of the cocoon (Equations 3 and 5)^. The normalization 
factors are based on the estimated ages, ambient gas densi- 
ties and 151-MHz flux densities for the FRII radio galaxies 
in the 7C Redshift survey. They also depend on the adopted 
valuc of a scaling factor / (as defincd in Willott et al. 1999) 
for the total energy inside the cocoon, inferred from the 
minimum-energy condition. The factor / parameterizes our 
ignorance of the exact geometry of the cocoon, low-energy 
cut-off in the electron distribution, filling factor, deviations 
from cquipartition, and proton contribution. The assump- 
tions we adopted in Section 3.6 to estimate the total energy 



* a value typicai of plasmas where electrons and positrons orig- 
inate as secondary particles foUowing coUisions of a primary 
hadronic jet with the ISM (Pacholczyk 1970). 



^ A similar method was followed by Heinz (2002) , see his Eq. (3) , 
but with the assumption of a radio spectral index «r = —0.5, 
definitely ruled out in our case. 
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correspond to a choice of / ~ 15, which is a plausible num- 
ber for QSO radio lobcs; it was estimated in Blundell & 
Rawlings (2000) that / is most likely to vary between ~ 10 
and 20. Therefore, our normalization coincides with the one 
adoptcd by Punsly (2005), that is £ = 4.0 x 10*^ crg for 
cPFi5i = 10^^ erg Hz~^ sr~^. This choice of normaliza- 
tion may easily vary by a factor of 2 at the reference flux, 
and laxger uncertainties may be introduced as we scaJe the 
solution from the parameter range suitable for radio-loud 
QSOs (c.g., n,H - 10"^ cm-'\ t - lO'^ yr, £. ~ 10*^ erg s'^), 
to the range suitable for microquasars (e. g., nn ~ 1 cm~^, 
t ~ 10^ yr, C < 10*^ erg s~^). Another caveat is that sources 
with a radio spcctral indcx ~ — 1 arc outlicrs in the FRll 
sample; most liavc indiccs between ~ —0.5 and ~ —0.9 (c.g., 
Figure 5 in Willott et al. 1999). 

Inserting the observed flux (0.28 mjy at 4.85 GHz) and 
speculative size (re « 30 pc) of the radio-emitting bubble 
into Equations (3-5), for an electron-positron plasma, we 
obtain an age t ^ 1.43 x 10® n]/^ yr, an average jet power £, = 
7.1 X lO^^njj^''^ erg s"\ a total energy Etot « 3.2 x 10®^ erg, 
and an expansion velocity v « (2rc/5t) « S2n^^^^ km s^^ at 
present time. We recall that we estimated (Section 3.6) a to- 
tal pressure Pin « 5x 10~^° dynes cm~^ inside the lobes; this 
is much larger than the thermal pressure of the undisturbed 
ISM, which is Pout ~ 3x 10"^^(riH/l cnr^)(r/10* K) dynes 
cm~^, with a sound speed ~ 15 km s^^ for T m 10"* K. 
The expansion of the bubble into the ISM is supersonic, 
with Mach number M w 82/15 « 5.5. Hence, we expect 
the expanding cocoon to drive a strong shock into tlic ISM. 
Pressure balancc should occur between the pressure inside 
the cocoon (Pin) and the thermal plus ram pressure of the 
shocked ambient gas, swept up in a shell outside the co- 
coon. Applying Standard Rankine-Hugoniot conditions (e. g., 
Spitzer 1978), we expect a thermal pressure in the shocked 
gas Ps = (5^4^ - l)Pout/4 « 37n-'Pout ^ 1 x 10"'° dynes 
cm~^, a density of the shocked gas ps w 3.6 times the am- 
bient ISM density, and a ram pressure « ps"^ w 4 x 10~^° 
dynes cm~^. Thus, even though our simple spherical model 
is only a crudc approximatioii of a real radio lobc, we ob- 
tain that the total pressure inside the cocoon is balanced by 
the thermal plus ram pressure of the swept-up interstellar 
medium, as we should expect. 

A higher jet power and shortcr bubble age are inferred 
if protons contribute a significant fraction of energy. For 
example, in the case of a 99% hadronic contribution, we 
obtain an age t w 5 x lO^'n^''^ yr, an average jet power 
C = 2 X lO^'^n^^''^ erg s~^, and an expansion velocity v w 
230ng^^^ km s~^ at present time. The total energies and 
pressure inside the bubble are an order of magnitude higher 
than in the leptonic case (--^ B^). The internal pressure is 
still roughly balanced by ram plus thermal pressure in the 
shocked gas (~ v^). 

Finally, from the Rankine-Hugoniot equations we also 
predict a temperature in the swept-out shell ~ 10 times the 
ambient ISM temperature, i. e. ~ 10'' K. The shell thickness 
should be « O.lrc ~ 3 pc. Wc infcr a cooling timescale for 
the shocked gas ~ 10*' yr, that is, louger than the estimated 
lifetimc of the radio cocoon. This may explain why we do not 
see stronger Balmer line emission from that region. The total 
synchrotron energy radiated over the lifetime of the bubble is 
also negligible (O(10~^)) compared with the injected energy 



Ct (Willott et al. 1999). This explains why the value of i?tot 
in Equation (4) can be treated as the total energy inside the 
bubble at present time. 

Considering the various sources of uncertainty men- 
tioncd earlier, the sccnario outlincd in this Section can 
only represent one possible solution in agreement with the 
(few) observed parameters. Nonetheless, it is a plausible 
scenario, consistent with existing models for microquasars. 
As already mentioned in Section 3.1, a steady radio jet is 
thought to dominate the power output in the LHS, and a 
transient jet may be present in the VHS, while it is sup- 
pressed in the HS. In the conventional definition of LHS, 
the X-ray luminosity Lx ^ 10~^LEdd, and the jet power 
Lj/LEdd ~ A(Lx/LEdd)'/^ with 6 X 10--' < ^ < 0.3 
(Fender, Gallo & .lonker 2003; Fcndcr et al. 2004; Malzac, 
Merloni & Fabian 2004). Hence, the jet power in the LHS 
is L J < 0.03LEdd. We compare this value with the value of 
jC ~ 7 X 10^* erg s~^ inferred from our self-similar soIutions. 
If most of the energy has been injected into the bubble by a 
steady LHS jet over the past ~ 10 ' yr, the Eddington lumi- 
nosity has to be > 2 X 10'*'' erg s~^, corresponding to a BH 
meiss > 150Mo. However, the upper limit to the luminos- 
ity in the LHS is somcwhat conventional: a steady jet may 
persist (with a similar scaling) as the source moves up in 
luminosity into the "hard" VHS, at Lx ~ L^dd (Fender et 
al. 2004) . If the source spent most of its time in that state, 
ifidd ~ a few 10^® erg s~* would sufRce. Finally, transient 
ejections with Lj ~ O.l-lLsdd occur at the spectral transi- 
tion between "hard" and "soft" VHS. 

We do not know how much time the X-ray source has 
spent in the different accretion states, and for what fraction 
of time it has been ofF. As we discussed in Section 3 (see 
also Soria et al. 2004), it is currcntly in somc kind of VHS, 
with Lx ~ lO**" erg s~^ and some ongoing Haring activ- 
ity suggesting coronal ejections. The estimated average jet 
power is at least consistent with a BH mass « IOOM0. For 
comparison, GRS 1915+105 {Mbh ~ MAf©) spends lorig 
periods of time in a flaring state with Lj ~ 10^® erg s~* 
and Lx w lO^'' erg s"* (Fender et al. 1999). The average jet 
power in SS433 is thought to be ~ 10^® erg s~^ (Fabrika 
et al. 2004, and referenccs therein). In the ULX class, a jet 
power ~ 1.5 x 10'^^ erg s^^ « 0.2Lx would be required to 
cxplaiu the ionized bubble around NGC 1313 X-2 (PakuU 
& Mirioni 2002). Another example — at the opposite end of 
the BH mass range — is the quasar PKS 0743—67, showing a 
radio jet with Lj ~ O.lLEdd ~ 10**^ erg s~^ as well as an 
accretion luminosity ~ L^dd (Punsly & Tingay 2005). 



5 CONCLUSIONS 

We have used the ATCA to study the radio counterpart of 
a bright ULX in the starburst irregular galaxy NGC 5408. 
The spectrum is rather steep (singlc power- law indcx ~ — 1), 
typical of optically-thin synchrotron. The flux has remained 
unchanged between March 2000 and December 2004. We 
speculate that the source is marginally resolved in a com- 
posite (4.8 -I- 6.2) GHz map, implying a size « 50 x 70 pc. 
However, deeper 6.2 or 8.4 GHz observations will be neces- 
sary to verify this hypothesis unequivocally, and to measure 
the spectral index more precisely. 

We have discussed the possibility that the source repre- 
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sents the coro radio omission from a ULX jet, as previously 
hypothcsizod (Kaarct ot al. 2003). We argued that a steady 
jet in the LHS is inconsistent with its steep radio spectrum, 
its very high flux (« 0.28 mjy at 4.8 GHz), and the X-ray 
spectral and timing proportics of the ULX. Thcrcforc, sim- 
ple BH mciss determinations based on fundamental-plane 
relations between LHS core radio and X-ray luminosities 
are not appUcable to this object. A sporadic, major radio 
outburst is also ruled out, given the constant radio flux be- 
tween the two cpochs of the ATCA obscrvations (2000 and 
2004), and the X-ray behaviour. Radio emission from a per- 
sistent, flaring VHS (perhaps similar to the dip-flare state 
in GRS 1915+105, characterizcd by a scqucnce of repcatcd, 
minor cjections) can be consistcnt with the observed flux, if 
enhanced by (moderate) relativistic Doppler boosting. We 
consider this scenario shghtly contrived, because there are 
currently no known examples of similar, persistent VHS flar- 
ing in Galactic microquasars. On the other hand, if the ac- 
creting BH has a mass ~ IOOM0 as suggcstcd by the X- 
ray luminosity (from the Eddington argumcnt), the spectral 
state transition pattern may be rather different from Galac- 
tic microquasars. The other reason why we consider core 
radio emission an unlikely scenario is that we have at least 
speculative evidence of a spatially resolved source. 

Instead, we suggest that the emission could be either 
from an underlying SNR (perhaps related to the formation 
of the compact object powering the ULX), or from radio 
lobes inflated by a relativistic jet from the ULX. A young 
shell-like radio SNR could reach sufHcient luminosity, would 
have a constant flux over many years, and has an optically- 
thin synchrotron spectrum. However, in contrast with its 
unusually high radio luminosity, this source is extremely 
faint in the optical. From a continuum-subtracted Ha im- 
age, one can see an H II region approximately coincident 
with the X-ray and radio sources. This region appears as 
a bright knot in a complex web of Ha-emitting filaments, 
connected to the main starburst complex. We estimate a 
luminosity ~ 1.5 x 10^^ erg s~^ in the line. This is an up- 
per limit to any possible contribution from an SNR. In fact, 
we showed that the few OB stars clearly detected in the 
same region, and/or the X-ray source itself would contributc 
enough ionizing photons to explain the Hq emission. A study 
of higher-ionization lines will be required to determine the 
ionization mechanism. Moreover, if we accept the evidence 
for a marginally-resolved radio source, Standard equiparti- 
tion arguments imply a total energy > a few x 10'''^ erg, and 
probably as high cis ~ 10^^ erg in the SNR: this is extremely 
high for a single SN, and may require a hypernova or a few 
normal SNe going off in the same group of stars. We con- 
clude that the SNR interpretation is unlikely, though not 
ruled out. 

A radio lobe powered by a jet is thus the most likely 
scenario, in our opinion. In equipartition, the total energy 
content is ^ 10®^ erg for a leptonic jet, or up to an order of 
magnitude higher if there is a significant hadronic contribu- 
tion. We modelled the lobe as a spherical bubble of radius 
30 pc, and solved a set of self-similar equations — often used 
in the literature to model radio lobes in FRII quasars — 
to relatc the bubble expansion with the energy input (i. e., 
the jet power) and the observed radio flux and spectral in- 
dex. Assuming an electron-positron plasma, we obtain that 
an average jet power w 7 x 10^* erg s~^ over a time-scale 



« 1.4x 10'' yr would indeed produee a 30-pc cocoon with the 
observed radio brightness, a total energy content ~ 3 x lO^'^ 
erg, with a randomly-oriented magnetic field «0.1 mGauss. 
Such a bubble expands into the ambient ISM with a present 
speed of ~ 80 km s~^, sweeping up a sliell of shocked gas. 
The pressure inside the bubble (« 5 x 10"^*^ dynes cm~^) 
is balanced by the thermal plus ram pressure of the gas in 
the swept-up shell. An age ~ 10® yr implies that the steep- 
ness of the spectrum around 5 GHz is from injection and 
cannot be duc to ageing (unless the magnetic field is well 
above equipartition); the age is also shorter than the radia- 
tive cooling timescale of the gas in the swept-up shell. 

An average jet power ~ 7 x 10^* erg s~^ is almost an 
order of magnitude higher than the jet power in the most X- 
ray luminous Galactic microquasar, GRS 1915-1-105, during 
its long-term VHS oscillation cycles. The corresponding X- 
ray source, NGC5408 X-l, is itself almost an order of mag- 
nitude more luminous than GRS 1915+105. On the other 
hand, an average jet power ~ 10'^'' erg s^^ has been inferred 
for SS 433 and for the ULX NGC 1313 X-2, both of which 
exhibit an extended lobe structure. We briefly summarized 
the accretion states in which a jet may be produced, and the 
relation between jet power, X-ray luminosity and Edding- 
ton luminosity. Our results are consistent with a luminosity 
Lx w I-Edd w 10*° erg s"^ (BH mass w lOOM©), and an 
average jet power Lj ~ 0.1Lx; they also suggest that the 
ULX has been at least moderately active (accretion power 
> 10~^LEdd) for most of the past 10® yr. 

Understanding this radio source in NGC 5408 is crucial 
for constraining the nature of ULXs, their association with 
starburst cnvironments, their process of accretion, and in 
particular the balance between radiative luminosity emit- 
ted by the accretion disk and mechanical luminosity car- 
ried the jet. Very few ULXs have a detected radio coun- 
terpart. Two ULXs in M 82 may have radio counterparts: 
one is probably a SNR, while the nature of the other one 
is still unclear (Kording et al. 2005); however, the physical 
association between those sources and the nearby ULXs is 
still uncertain, given the density of sources in that region. 
The only other unambiguous ULX-radio source association 
is in Holmberg II (Miller, Mushotzky & Neff 2005). In that 
case, too, the spectrum is consistent with optically-thin syn- 
chrotron, though less steep than in our case {an ^ —0.5). 
The Holmberg II radio source is resolved (size w 50 pc) and 
has a flux density « 0.7 mJy at 5 GHz (Miller et al. 2005), 
corresponding to ~ 0.3 mJy at the distance of NGC 5408; 
hence, the two sources may be very similax (the X-ray lu- 
minosity and spectral properties of the corresponding ULXs 
are also similar: Dewangan et al. 2004; Soria et al. 2004). 
In that case, the purc SNR scenario appears inconsistent 
with the optical line spectrum; the most likely interpreta- 
tion is that the synchrotron bubble is powered by the ULX 
jet (Miller et al. 2005). 

Finally, should we expect to find many more micro- 
quasar lobes in nearby galaxies, perhaps misidentificd as 
SNR? Taking NGC 5408 and Holmberg II as typical cases, 
such lobes would require a jet injection power > 10'^'' erg 
s~^ over > 10'' yr, in order to be detected by existiiig radio 
surveys of nearby galaxies. If we assume that the jet power 
L J ~ Z/x < ^Edd, such Systems should be associated with 
potential ULXs. In some cases, the X-ray source may not 
be currently active (or active at ULX level), yet there would 
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still be a bright synchrotron radio source, formed by the in- 
tegrated source activity over its lifetime. Therefore, finding 
the number ratio between microquasar lobes with or without 
a currently active ULX core could in principle provide a con- 
straint on the average duty cycle of ULX microquasars. As 
an example, two bright radio sources, consistent with syn- 
chrotron bubbles with minimum-energies > 10^^ erg (hence, 
candidate microquasar lobes or hypernovae) have recently 
been discovered in NGC 7424 (Soria et al. 2006) . One is as- 
sociated with a transient ULX, the other has no correspond- 
ing X-ray source, down to a completeness limit ~ 3 x lO^'^ 
erg s~^. Thus, one could speculate that the ULX responsible 
for latter radio source is currently in the off state. Further 
discussion of these issues is beyond the scope of this paper. 
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APPENDIX A: NOTES ON THE 
MINIMUM-ENERGY CONDITION 

Two alternative formulations are used in the literature to es- 
timate the minimum energy associated with a radio source 
(e.g., Tudose et al. 2006). One choice (e.g., Pohl 1993, Bick- 
nell 2005), arguably the more physical, is to assume an en- 
ergy range (7min, 7max) for the relativistic electrons. In this 
framework, following Bicknell (2005), a useful way of ex- 
pressing the minimum- energy magnetic field is: 
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where the energy spectrum of the electrons is N{E)dE ~ 
E~^dE, or = (1 — p)/2, me and e are the electron mass and 
charge, c the speed of light, the specific surface bright- 
ness (flux per angular beam area), (2rc) is the characteristic 
source diameter, and k is the ratio between the energy of 
protons and electrons. We have assumed a filling factor of 
1, for simplicity. The functions 
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and r(z) is the Gamma function. The corresponding total 
(minimum) energy density is (Bicknell 2005): 
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In our analysis, we adopted 7inin ~ 1, 7max ~ 10*. For 
a steep spectrum, the minimum energy depends only very 
weakly on the high-energy cut-off. As for the low-energy 
cut-off, there are theoretical arguments and observational 
evidence for the presence of low-energy electrons in a radio 
lobe (e.g., BlundeU & Rawlings 2000, their Section 11.1). 
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The alternative choice (e. g., Burbidge 1956; Pacholczyk 
1970; Longair 1994; Willott et al. 1999) is to assume a fre- 
) of the radio emission. In that case, 
the magnetic field for which the energy density is minimized 
is (Longair 1994): 



« 119.7 C2(p)*/^(l + fe)^/^(^) 
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where V is the volume of the bubble, and we have assumed 
a unit filling factor. The functions 
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and T{z) is again the Gamma function. The minimum en- 
ergy density is then (e. g., Tudose et al. 2006): 
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The physical connection between the two formahsms is 
that for an electron of energy 7, the peak of its synchrotron 
emission is at a frequency 
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In some cascs (c. g., Willott et al. 1999; Punsly 2005), the 
minimum-energy quantities arc ostimatcd for a conventional 
lower cut-off at i^min = 10 MHz (lowest frequency accessi- 
ble to Earth-based radio observations), and then the energy 
density is corrected by a numcrical factor / > 1 to account 
for low-encrgy olectrons emitting below 10 MHz, and for 
hadronic contributions. In our case, setting a low-frequcncy 
limit of 10 MHz would imply an implausibly high 7min ~ 300 
for a leptonic plasma. Empirically, it was shown (Blundell 
& Rawhngs 2000; Punsly 2005) that the energy in the lobcs 
of radio galaxics is ~ 10-20 times the minimum- energy esti- 
mated from a leptonic plasma with a 10-MHz cut-off. For a 
radio spectral index w —1, our choice of a leptonic plasma 
with 7min ~ 1, 7max 10* corresponds to a correction factor 
/ w 15, consistent with the observations (see also Section 4) 



where the density p(r) ~ por^'' . For a constant ISM density 
(( = 0), Eq.(Bl) reduces to Eq.(3). The bubble size is also 
related to the total energy Etot = Vetot ~ r^etot, and from 
Eq.(4): 



In the minimum-energy approximation, etot = tn 
from Equations (A1,A5), 

With simplo algcbraic manipulations, the system of equa- 
tions (B1,B2,B3) can be solved as: 
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which is identical to Eq.(5) for ^ = O and p = 3. Similarly to 
Heinz (2002) (who discussed the case p = 2), the normaliza- 

tion coefScicnt in front of Eq.(B4) was fixed to be consistent 
with Willott et al. (1999) 's study of radio lobe properties in 
radio galaxies (Section 4) . 



APPENDIX B: TIME EVOLUTION OF THE 
RADIO BUBBLE 

Wo outhne here the derivation of Equation (5) from Equa- 
tions (3,4) and the minimum-energy condition. Adopting 
Falle (1991)'s self-similar solution, the radius Vc of an ex- 
panding cocoon scales as 
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